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Abstract 

We investigate static space dependent cr(x) = (•ipip) saddle point con- 
figurations in the two dimensional Gross-Neveu model in the large N limit. We 
solve the saddle point condition for <j(x) explicitly by employing supersymmet- 
ric quantum mechanics and using simple properties of the diagonal resolvent of 
one dimensional Schrodinger operators rather than inverse scattering techniques. 
The resulting solutions in the sector of unbroken supersymmetry are the Callan- 
Coleman-Gross-Zee kink configurations. We thus provide a direct and clean 
construction of these kinks. In the sector of broken supersymmetry we derive the 
DHN saddle point configurations. Our method of finding such non-trivial static 
configurations may be applied also in other two dimensional field theories. 
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I Introduction 



The Gross-Neveu model JT| is a well known two dimensional field theory of iV massless 
Dirac fermions tp a (a = 1, . . . , N) with U(N) invariant self interactions, whose action 
is given by 



S = d 2 x 



N g 2 / N y 

a=l A \a=l / 



We are interested in the large N limit of (|I|) in which N — > oo while Ng 2 is held fixed. 
Decomposing each Dirac spinor into two Majorana spinors one observes that 5* is 
invariant under 0(2N) flavour symmetry containing the U(N) mentioned above as a 
subgroup [[J. The field theory defined by ([!]) is a renormalisable field theory exhibiting 
asymptotic freedom, dynamical symmetry breaking and dimensional transmutation. 
Its spectrum was calculated semiclassically (in the large N limit) in 0. It contains 
the fermions in ([]]) which become massive, as well as a rich collection of bound 
states thereof (the so called "DHN states"). The spectrum of (|l|) contains also kink 
configurations 0, [5]. We refer to these as the Callan-Coleman-Gross-Zee (CCGZ) 
kinks in the sequel. These kinks are expected to be part of the spectrum of the 
Gross-Neveu model since dynamical breaking of the discrete chiral symmetry in the 
Gross-Neveu model suggests that there should be extremal field configurations that 
interpolate between the two minima of the effective potential associated with (|T|) in 
much the same way that such configurations arise in classical field theories whose 
potential term has two or more equivalent minima. 

The Gross-Neveu model has also a system of infinitely many (non-local) conserva- 
tion laws which (presumably) forbid particle production in scattering processes and 
enables the exact calculation of S"-matrix elements in the various sectors of the model 
. Results of such calculations are in agreement with the "Large N" calculation of 
the spectrum. 

In this paper we discuss static space dependent a(x) = (ipip) configurations that 
are solutions of the saddle point equation governing the effective action corresponding 
to ([[]) as iV — > oo. Such a(x) configurations correspond to non-trivial excitations of 
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the vacuum [7], |8|] and are therefore important in determining the entire spectrum of 
the field theory in question^] and its finite temperature behaviour as wellp. Such 
configurations are important also in discussing the behaviour of the expansion of 
© at large orders §0,g. 

Our discussion makes use of supersymmetric quantum mechanics and simple prop- 
erties of the diagonal resolvent of one dimensional Schrodinger operators. Using these 
two basic tools we are able to solve the saddle point condition for static cr(x) config- 
urations explicitly. The supersymmetry alluded to above relates the upper and lower 
components of spinors, implying that the square of the Dirac operator may be decom- 
posed into two isospectral Schrodinger operators. It is closely related to the soliton 
degeneracy discussed in where the soliton is considered as a degenerate doublet 
having fermion numbers ±|. Our explicit solution of the saddle point condition in 
the sector of unbroken supersymmetry consists of the CCGZ kink configurations and 
it therefore provides a clean and direct construction of these kinks. In the sector 
of broken supersymmetry, our explicit solution reproduces the so called DHN saddle 
point configurations p[ . 



In a recent paper [13], we have applied a similar method to the anharmonic O(N) 
oscillator and the two dimensional O(N) vector model in the limit iV — > oo. In the 
latter case we have found that the effective action sustains in the large temperature 
limit extremal bilinear condensates of the O(N) vector field that are analogous to the 
CCGZ kinks in the Gross-Neveu model. 

The paper is organised as follows: In Section II we analyse the saddle point equa- 
tion for static <j{x) configurations employing supersymmetric quantum mechanics. 
Using simple manipulations, we show that the latter equation may be expressed in 
terms of the hamiltonian of only one of the supersymmetric sectors. In Section III. 1 
we resolve the saddle point equation into frequencies and demand that the extremum 
condition be satisfied by each Fourier component separately. This strong condition 
turns out to leave us always in the sector of unbroken supersymmetry. This corre- 
sponds physically to cr(x) configurations that interpolate between the two vacua of 
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the Gross-Neveu model at the two ends of the world, which are the CCGZ kinks that 
we find as explicit solutions of the saddle point equation. In Section III. 2 we solve the 
saddle point equation without separating it into frequencies, under the assumption 
that there is only one bound state. This leads to the DHN a(x) configurations which 
belong to the sector of broken supersymmetry. We conclude our discussion in section 
IV. 
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II The Saddle Point Equation for Static Solutions 



Following]]]] we rewrite ([!]) as 

S = d 2 x 



gaipip a z 



(2) 



where o~(x) is an auxiliary field [ p~5|j . 

Thus, the partition function associated with (0) is 



Z 



T>a Vip Dip exp < i 



oo oo 

dt / dx 



tjji0ip — gcnptp — - a 2 



(3) 



Gaussian integration over the grassmannian variables is straightforward, leading to 
Z = J T>o~ exp{iS e ff[a}} where the bare effective action is |T6[] 



S effW] 



dt I dxa 2 — i — Tr In 



ga)(i0+ gar) 



(4) 



The ground state of the Gross-Neveu model (|l|) is described by the simplest extremum 
of S'e/zl] in which a = cr is a constant that is fixed by the (bare) gap equation 

d 2 k 1 



-go- + iNg 2 tr J 



(2tt) 2 #-g(T 



0. 



(5) 



Therefore, the dynamically generated mass of small fluctuations of the Dirac fields 
around this vacuum is 

l— " 

3 Ng 2 „M 



m = ga = lie 



(6) 



where [i is an arbitrary renormalisation scale, and the renormalised coupling gn (//) 

* l J — 

is related to the cut-off dependent bare coupling via A e N a 2 ( A ) = fie N3 * M where 
A is an ultraviolet cutoff. Since m is the physical mass of the fermions it must be 
a renormalisation group invariant, and this fixes the scale dependence of the renor- 
malised coupling gn ,namely, equation @. From now on we will drop the subscript 
R from the renormalised coupling and simply denote it by Ng 2 . 

As was explained in the introduction, we are interested in more complicated ex- 
trema of S e ff, namely static space dependent solutions of the extremum condition on 



S e ff. This condition reads generally 

SS eff 



5<t(x, t) 
N 



-<j(x, t) 



- t -tr{[2g 2 a{x,t)+tgYd^ (x, t\ [ □ + g 2 a 2 - igj^a]' 1 ^, t)} = (7) 

where "tr" is a trace over Dirac indices. 

Specialising to static cr(x) configurations, and using the Majorana representation 
7 1 = ia 3 , 7 = cr 2 for 7 matrices, (0) becomes 



2i 



-a(x) = tr 



2ga — d x 
2ga + d x 



du ( R + (x,lu 2 ) 
2tt{ R-{x,u 2 ) 



dui 
2^ 



(2ga - d x ) R + (x, to 2 ) + (2ga + d x ) f2_ (. 



x, UJ 



(8) 



where 



R±(x, uj 1 



h+ — co> 2 



are the diagonal resolvents of the one dimensional Schrodinger operators 



h ± = -d 2 x +g 2 a 2 ±ga\x) 



(9) 



(10) 



evaluated at spectral parameter uj 2 



Note that h± are positive semidefinite isopectral (up to zero-modes) hamiltonians, 
since (|TUp may be brought into the form |T9], M 

/i + = Q^Q , /i- = QQ^ where 

Q = ~ + go , Q f = + per. (11) 

These operators may be composed into a supersymmetric hamiltonian H = 
describing one bosonic and one fermionic degrees of freedom Hl9| , which we identify 
with the upper and lower components of the spinors. Supersymmetry implies here 



( h + \ 
(0 hi 
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that an interchange of the bosonic and fermionic sectors of H leaves dynamics un- 
changed, as can be seen from the fact that Eqs. (|^)-(|ITD are invariant under the 
simultaneous interchanges 



a — > —a , h± — > h T , 

R± -> R T . (12) 

Isospectrality of h + and h- alluded to above means that to each eigenvector ip n 
of h- with a positive eigenvalue E n , there is a corresponding eigenvector n of h + 
with the same eigenvalue and norm, and vice-versa. The precise form of this pairing 
relation is 

1 ~t. 



ij n = -^=Q(j) n ; E n >0. (13) 



It is clear that the pairing in pB] ) fails when E n = 0. Thus, in general one cannot 
relate the eigenvectors with zero-eigenvalue (i.e.-the normalisable zero-modes) of one 
hamiltonian in ([TOD to these of the other. Should such a normalisable zero-mode 
appear in the spectrum of one of the positive semidefmite operators in (|10D, it must 
be the ground state of that hamiltonian. In this case the lowest eigenvalue of the 
supersymmetric hamiltonian H is zero, which is the case of unbroken supersymmetry. 
If such a normalisable zero-mode does not appear in the spectrum, all eigenvalues 
of H, and in particular-its ground state energy, are positive, and supersymmetry is 
broken. Since the ground state of a Schrodinger operator is non-degenerated, h± can 
have each no more than one such a normalisable zero-mode. Moreover, it is clear 
from (III]) that in our one dimensional case, only one operator in (|T0|) may have a 
normalisable zero-mode, since it must be annihilated either by Q or by Q\ In cases 
of unbroken supersymmetry, we will take such a normalisable zero-mode to be an 
eigenstate of h-, namely, the real function 

x 

$ (x) = Me o (14) 



which is the normalisable solution of the differential equation 

g f ^o = o (15) 

where M is a normalisation coefficient. 

Note that a necessary condition for the normalisability of \l/o is that <j(x) have 
the opposite behaviour at ±00. Thus, physically, cases of unbroken supersymmetry 
lead to cr(x) configurations that interpolate between the two vacua of (J4|) at the two 
ends of the world, while cases of broken supersymmetry yield cr(x) configurations that 
leave and return to the same vacuum state. 

Assigning the zero-mode to h- poses no loss of generality, since the other possible 
case is related to this one via (|T2| ) . In what follows we will denote the right hand side 
of fllil) by $0 a l so m cases of broken supersymmetry where it is non- normalisable. 
This should not cause any confusion, since the ground state will be denoted by if) , 
which will be equal to when supersymmetry is unbroken. 

By definition, the diagonal resolvents in (||) are given by the eigenfunction expan- 
sions 

R+ {x) = g v^ML (16) 

n>0 - U 

where the sums extend over all eigenstates, including the continua of scattering states 
where they are understood as integrals. 

Using ( [T3D , R + may be expressed in terms of the ^'s as 

R + (X,U ) = £ - ^j-- - 



2 

„ E n ml{x) E n -u* (17) 



where W n (x) = ^ (x) , ijj' n (x) — ^' (x)ijj n (x) is the wronskian of and ip n . An ele- 
mentary consequence of the Schrodinger equation is that W n (x) = —E n ty (x)ip n (x). 
Using this relation, (|16D and (|17]) imply the important relation 



( 2,,a c) r ) R + - (2ga + d x )R_ = 2[ 2ga + A j ( x | 7 ^ i_ \ x ) . 



Here V is the projector 

V = l-A|^o)<V>o| (19) 

that projects out the ground state of 7t_ when supersymmetry is unbroken (A = 1), 
and is just the unit operator otherwise (A = 0). We can also use (|5|) to make a 
frequency resolution of unity as 

A 



Ng* 



duo 



71 



(x\ 



1 



-&i + m? — uj 2 — ie 



\x) . 



(20) 



Moreover, from (|T8|), (|20| ) and the elementary relation 

1 , , 1 



(x\ 



—d 2 + m 2 — uj 2 — ie 



\x) 



2\Jm 2 — uj 2 — ie 



(21) 



the frequency resolution of (|]) becomes 
1 



Tx^ v 



h_ — ur 



\x) = 2ga(x) 



1 



- (x\V 



h_ — 



\x) 



1\/m 2 — uj 2 

Substituting ( |19|) into the last equation, all dependence on A cancels out and the 
extremum condition (Rf) is shaped into its final form 



du j dR^ (x, uj 2 ) 
2^ 1 dx 



2ga(x) 



2v^ 



(22) 



which has to be satisfied regardless of whether supersymmetry is broken or not. 
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Ill Static Solutions to the Extremum Condition 
IIL1 Case of Unbroken Supersymmetry 

The simplest way to look for solutions of the static saddle point equation (^) is to 
demand that it be satisfied by each frequency mode separately, namely restricting i?_ 
by the differential condition 



dR_ (x, k? 2 
dx 



2ga(x) 



— R- (x, uj 



(23) 



Now, i?_ (x, uj 2 ) , being the diagonal resolvent of h_ at spectral parameter uj 2 is 



subjected to the so-called " Gel'fand-Dikii" equation |p^|||18| 



- 2R_(x, uj 2 )R'^(x, uj 2 ) + [R'_(x, uj 2 )) + AR 2 _(x, uj 



2 2 

9 ° 



go' - uj 2 



(24) 



Therefore, both equations ( p3[ ) and ( Ell ) must hold and they form a system of 
coupled non-linear differential equations in the unknowns <j(x), and i?_(x,aj 2 ). Sub- 
stituting R'_ and R"_ from (|23|) into (p4]) we obtain a quadratic equation for R_ whose 
solutions are 



R^ (x, uj 2 ) 



-go' ± J(go'Y + Au 2 g 2 o 2 - Alu 2 (m 2 - uj 2 ) 



Auj 2 V' 



(25) 



To see what the two signs of the square root correspond to we observe that the 
solution with the negative sign in front of the square root has a simple pole as a 
function of uj 2 at uj 2 = with a negative residue, while the other solution is regular 
and positive at u 2 = 0. Therefore, from fliTf ) it is clear that the negative sign root 
corresponds to the case of unbroken supersymmetry, where the simple pole signals the 
existence of a normalisable zero-mode in the spectrum of h-, while the positive root 
solution corresponds, for similar reasons, to cases in which h_ lacks such a zero-mode. 
We will see below that the latter solution corresponds also to the case of unbroken 
supersymmetry, where the zero-mode is in the spectrum of h + . 



We concentrate now on the branch cut singularity in (|25|) at uj 



2 which 



signals the threshold of the continuous part of the spectrum of h_. Clearly, h_ can 
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have no continuous spectrum other than that starting at uj 2 = m 2 , thus leading to 
the conclusion that _R_ in fl25|) can have no other branch points in the uj 2 plane. 
Therefore, the expression under the square root in (^) must be a perfect square as 
a polynomial in u 2 , namely, 

± ga \x) = g 2 a 2 (x) - m 2 . (26) 

From (p6|) we find straight forwardly the solutions 

ga(x) = ± mtanh [m (x — x )} (27) 

which are exactly the CCGZ kinks and anti-kinks. Here the parameter xq is an inte- 
gration constant that implies translational invariance of ([27]) since it is the arbitrary 
location of the kinks. Clearly, both cases in (^), and therefore both cases in 
lead to h± operators that do not break supersymmetry, since (|14]),(^]) imply that 



i ,.i 



^ = 2 e -9j x °(y)dy = 2 S ech[m(x - x Q )] (28) 

is the normalisable zero-mode of /i_ for the kink configuration, and of h + when a(x) 
is the anti-kink. Due to zero binding energy, fermions trapped in this potential do 
not react back on the u{x) field [0. 

Note that in deriving ([27]) above we have not set any apriori restrictions on i?„, 
thus, it is a very interesting question whether ( p7|) are the only possible static extremal 
u{x) configurations in the sector of unbroken supersymmetry or not. 



As it stands, our frequency decomposition of the extremum condition (23) seems 
to lead always to extremal u{x) configurations that do not break supersymmetry. We 
are thus unable to find in this manner the extremal cr{x) configurations found in 
in which o{x) has the shape of kink anti-kink pair that are very close to each other. 
Such a configuration evidently breaks supersymmetry Jl9[ . They will be discussed in 
the next subsection. 

When cr(x) is a kink (|25|) becomes 

R.^u 2 ) = - m " ^ 2 M^o)] + 1 (29) 
2u; 2 vm 2 — uj 2 2ym 2 — uj 2 
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while for anti-kinks it is just the expression on the right hand side of 

These statements on R_ are consistent with the explicit form of the hamiltonians 
h±. Using ([TTJI) , (P7|) we find in the kink case 

h + = -d\ + g 2 o 2 + go' = -d 2 x + m 2 and 

h- = —d 2 + g 2 a 2 — ga' = —d 2 + m 2 — 2m 2 sech 2 [m(x — x )} (30) 

while in the anti-kink case h± interchange their roles . Thus, in the latter case, /i_ 
becomes the Schrodinger operator of a freely moving particle in a constant potential 
m 2 which is the reason why i?_ is given by the simple expression in fl2l|) in the anti- 
kink case. Indeed, that expression is the x independent solution of the Gel'fand-Dikii 
equation (|23|) corresponding to the constant potential g 2 o 2 — go' = m 2 that is positive 
at J 2 = 0. 

In the kink case, we have obtained the explicit form fl29f ) for i?_ by substituting 
( |27| ) into (p5[). As an independent check, we may deduce this expression for the 



diagonal resolvent for the potential g 2 o 2 — go' = m 2 — 2m 2 sech[m(x — xo)] in other 
ways. The simplest one is to apply an ansatz of the form i?_ = a sech 2 (/3x) + 7 to 
the Gel'fand-Dikii equation. Another way to derive it is to use the general formula 
R(x) = G(x,y)\ x=y = ^w^lp^ \ x=y where the ^'s are the so called Jost functions 
of the problem. In this case they are hypergeometric functions multiplied by sech 
factors. The wronskian in this expression is a ratio of T functions dependent on uj 2 
and m 2 [g]. 

We can deduce from (|25D the CCGZ solution (|27|) to the extremum condition 
(|25|) , (|2^) by yet another method which does not rely upon the brach-cut structure of 
i?_ but rather on its pole structure. Considering the case where the zero-mode is in 
the spectrum of h- we make a Laurent expansion of the appropriate expression for 



i?_ in (^5]) around the simple pole at uj 2 = 0. The leading term in this expansion is 

(31) 



1 

R- ~ 

UJ 1 



9 -^ + oU 

2m V 
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Comparing ([H]) to ([H]) and fllB]) we find 

^(a;) = M 2 ex P (-2g J a(y)dy\ = |£ (32) 



which yields 



ui 



x 



AfV 2 * = ^ ; <P(x) = g I a(y)dy . (33) 



o 



Solving (|33"D we find 

<j>(x) = lncosh[m(x — xq)\ + c (34) 
where c and xo are integration constants and we have imposed the normalisation 

oo 

condition / ^ 2 = 1 to fix W = Jf c . Clearly, then, f|) yields the CCGZ kinks 
upon differentiation. 

We are now in a position to verify briefly that the CCGZ kink configuration leads 
to an h- operator that has indeed a single normalisable zero-mode as the explicit 
form fl29|) for R_ suggests. Our discussion follows p0| . In the kink sector (|30| ) implies 



that the eigenstates of h + are simply these of freely moving particles (f>k{x) = e* fcx , 
with a continuum of strictly positive eigenvalues E + = k 2 + m 2 > m 2 . These states 
are isospectral to the eigenstates of h_, 

■ipk{x) = Q4>k{x) = { { e (35) 

which are therefore the scattering states of The S matrix associated with h_ is 
thus 

7 h — rn / h \ 

(36) 



ik-m 

bik) = — = exp 

ik + m 



. / k" 
i 7T — 2 arctan — 



Since /i + in (^) has only scattering states, 7t_ can have no bound states other than 
its zero-mode fl2"8p which must therefore be its ground state. This single normalisable 
state of corresponds to the single pole of S(k) in (|36|) at k — im. Note further 
that there are no reflected waves in any of the scattering eigenstates (|35|) of the 
Schrodinger operators in (|30D . This is also the case for the supersymmetry breaking 
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a(x) configurations in [Q] as well as in other exactly soluble models in two space-time 
dimensions [0. 

The fact that h± evaluated at the extremal point must be reflectionless can be 
deduced even without solving the extremum condition (^) explicitly, provided one 
makes apriori an assumption that /i_ has only a single bound state (namely, its ground 
state) regardless of whether supersymmetry is broken or not. 

To this end we consider (|23|) , in which i?_ may be replaced by Rp = (x\ V , _ 2 |x) 
as mentioned in the discussion preceding fl23|). Since h- is assumed to have a single 
bound state, Rp contains only scattering states of whose corresponding continu- 
ous eigenvalues E = k 2 + m 2 start at E = m 2 . We may therefore write the spectral 
resolution of Rp as 

flp= L4£MI! = L/c »W , (37) 
J E — ur J k z + m z — or 

where (x) = 27r\k\\ipk (x) \ 2 . In terms of pk the extremum condition ( [23] ) becomes 

d_Pk{x)_ = d 1 
rfx (x) <ix (x) 

whose general solution is 

p fe (x) = 1 + c k m 2 (x) (39) 

where is an integration constant. Since by definition pj, (x) cannot blow up at 
infinity, we must set = in the case of broken supersymmetry. Whether su- 
persymmetry is broken or not p k (x) obviously obtains the asymptotic value of 1 as 
x — > ±oo. Therefore, the scattering states ipk {%) of h_ are given by 

i) k (x) = yfl + c k V 2 (x) e ia ^ (40) 

where a k is a real phase. Substituting these functions into the eigenvalue equation for 
/i_ and considering its asymptotic behaviour as x — > ±oo we see, using the boundary 
condition go (±oo) = ±m that the phase becomes that of a free particle, which is 
obvious, but unimodularity of the phase factor implies further that there be only right 
moving or only left moving wave in ip k (x). Therefore, h_ must be reflectionless. 
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The physical significance of the CCGZ kinks is as follows: As was mentioned in 
the introduction, the dynamical properties of (H) are consequences of the fact that the 
("large N") effective potential V e g(a) extracted from (SI) has two symmetric equivalent 
minima at (er) V ac = ±o"o 7^ 0. This causes a dynamical breakdown of the discrete 
(Z2) chiral symmetry of (fj) under the transformation ip — > 75^, a — > —a , where 
the fermions fluctuating near the (a") vac = ±<To ground state acquire dynamical mass 
m = iperoU- In a similar manner to the appearance of kinks in classical field theories 
with potentials exhibiting spontaneous symmetry breaking, one should expect similar 
configurations to appear in field theories whose effective potential implies dynamical 
symmetry breaking. The CCGZ kinks (anti-kinks) are precisely such static space- 
dependent <j(x) configurations that interpolate between the two minima of V e g(a), 
and our calculations provide an explicit proof that they are indeed extremal points 
of (3). The various states appearing in the background of the CCGZ kink may be 
deduced by calculating the "partition function" of the Dirac field fluctuations in 
that specific a(x) configuration for a finite time lapse T (i.e.-the trace over the time 
evolution operator e~ lHT ) . This has been done explicitly in 0. The result is 

. T oo 

2~ l ^ i ^ X ( C 

-00 

exp jiiVT [<J k ink] - Yl Ui [a ]j - inu b [a kinh ] T j . (41) 

In this equation n is the total number of fermions and anti-fermions that are trapped 
in the single bound state of the kink, Ui [a] = \TE~i is the energy of the i —th state 



^ —iHT ^ (2 AO' 

Tre = > ——-7 rrexp 



■ °kink a 



in the background of a[|2l] and Ub is the energy of the bound state, which is zero for 
the CCGZ kinks. Using this and also the fact that in this background h± in ( p0|) are 
isospectral, we see that all terms in the second exponent in ([y]) are cancelled, leaving 



only the first exponent which is the mass of the kink, namely p2|, 

Mkink = z - I 17 n - <J 



2 



kink 



mN , , 

(42) 



Therefore, we see that all states contributing to (|4T|) are degenerate in energy (all 
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having the kink mass as energy) which is a direct result of the fact that out, = 0. Clearly 
there are 2 2N states in all, that form a huge reducible supermultiplet of 0(2N). Its 
decomposition into irreducible components is clear from the combinatorial prefactors 



in (f|l]) that simply tell us that the various bound states in the kink fall into anti- 
symmetric tensor representations of 0(2N) (where the integer n is the rank of the 
tensor) [0]. 

We close this subsection by checking explicitly that the CCGZ kink configurations 
obtained above indeed extremise the effective action in (^). Substituting the kink 
configuration in fl27f) and the explicit expressions (|29|) and (|2l|) for R_ and R + into 
the extremum condition (||) we find that the pole at u 2 = disappears from the right 
hand side of (§) in accordance with ([18]) leaving in the sum over frequencies only 
contributions from the scattering states. Thus, Eq.(|8|) becomes 



1 + iNg 2 J 



duo 



2n y/m 2 — uj 2 — ie 



ga(x) 



0. 



(43) 



implying that the term in the square brackets on the left hand side must vanish. But 
vanishing of the latter is precisely the Minkowsky space gap-equation of the Gross- 
Neveu model (||) for the dynamical mass m and it must therefore hold, confirming 
that the kinks in fl27T) are indeed solutions of the extremum condition (|8|). 
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III. 2 Case of Broken Supersymmetry 

Enforcing the saddle point condition at each frequency component in (|23|) led us 
directly to the sector of unbroken supersymmetry without any further assumptions 
on u) 2 ). Therefore, in order to find static extremal cr(x) configurations that lead 



to supersymmetry breaking, we must solve (22) as a whole. At a first sight this seems 



to be unmanageable p3| , since we apparently cannot use the Gel'fand-Dikii equation 
fl24D for i2_ which was so crucial for our treatment in the previous section. However, 
assuming (as in ||) that the required <j(x) configuration yields an h- operator with 
a single bound state at positive energy E b = uj 2 < m 2 , in addition to the obvious 
continuum of unbound fermions of mass m, our experience gained in the previous 
section leads us to the most general form of R_(x,uj 2 ) consistent with (H) and the 
preassumed form of the spectrum of This generic form of i?_ contains enough 
information in order to solve ([2"2"|). 

In order to construct this generic form of we note the following points: 

1. R- has scale dimension —1 in mass units as can be seen from its definition (|9]) 
or from the explicit expressions (E5f) and fl29|). 



2. Clearly, gu{x) must reach asymptotically either one of the vacua, namely, 
|gcr(±oo)| = m, ga'(±oo) = (actually, since we will end up indeed in the 
sector of broken supersymmetry, the two asymptotic values of go will turn out 
to be equal). Thus, the h_ operator resulting from such a ct{x) configuration 
must have the asymptotic form —d 2 + m 2 as \x\ — > oo. Correspondingly, its 



resolvent must have the asymptotic form ([2l|). 



3. The assumption that h_ has a single bound state at positive energy < uj\ < m 2 
implies that the analytic structure of R_(x,u> 2 ) in the complex uj plane must 
include in addition to the branch cuts along the real rays (— oo, —m) and (m, oo) 
only two simple poles on the real axis at uj = ±a>b . 

4. Any ansatz for R_ must obey the Gel'fand-Dikii equation (pjj). 



16 



Points number 1 and number 3 lead immediately to the general form 

(a?5 — uj 2 ) \/m 2 — uj 2 

where A,B,C and D are dimensionless numbers yet to be determined. Point number 
2 imply then that C = ~ and S = — £). 

In order to comply with point number 4, we substitute flUj) into the Gel'fand-Dikii 
equation (p4|). Multiplying the resultant expression through by (u 2 —u 2 ) 2 (m 2 —uj 2 ) we 
obtain an even quartic polynomial in cj|24| which must vanish. This yields three non- 
linear differential equations for go~(x) stemming from nullification of the coefficients 
of cu 4 , uj 2 and uj° in that polynomial. The condition which results from setting the 
coefficient of u 4 to zero reads 

(4A - l)(ga)> + (1 - AD)[{ga) 2 - m 2 ] = (45) 

which can be used to eliminate A and D from (f44|). Doing so (44) becomes 



R-(x,u ) = —— ^ + (46) 

4 (a;^ — uj z ) y/m z — u; z l\Jm A — co> z 

This expression is evidently very similar to expressions (|25|) and (|29D , the only differ- 
ence being that the double pole at uj = in those equations is resolved here into the 
two simple poles at uj — ±uv 

We must further subject ( |4"E| ) to the two remaining differential equations men- 
tioned above. It turns out that the equation resulting from setting the coefficient 
of u 2 to zero is simply the derivative (with respect to x) of the equation associated 
with the coefficient of a; [25 1. The latter is the differential equation for gu{x) we are 



looking for, which must be solved subjected to the boundary conditions mentioned 
in point number 2 above. Equivalently, and this is what we do below, we can read off 
the expression for the normalised bound state wave function i^b(x) in terms of gcr{x) 
from (f4"6]). The Schrodinger equation for ipb{%) provides than the required condition 
for gcr(x). 
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As can be seen from ([16]) and in a similar manner to our discussion following (|3T| ) 
in the previous subsection, we identify the residue of the simple pole of R(x, E = u 2 ) 
in (M) at E = u 2 as -^{x). Thus, 



go 1 — (gcr) + m 2 



^ 4yW -iu 2 ^4 



m 2 — V 



m 2 — cu? 



where 



V(ar) = g 2 a 2 - ga' 



is the potential of h-. Imposing the Schrodinger equation on (|47|) we have 



d x + V ( X ) \/ m2 - V ( X ) = ^V m2 - V ( X ) ■ 



(47) 



(48) 



(49) 



The solution to this equation, compatible with boundary conditions at infinity is 



V(x) = m 2 — 2K 2 sech 2 [n(x — Xq)} 



(50) 



where k 2 = m 2 — io\ and Xq is an integration constant. The corresponding bound 
state and resolvent are therefore 



-sech[n(x — x )] 

K 2 Sech 2 [K(x — Xq) 



(51) 



2 (a; 2 — c<j 2 ) \/m 2 — u 2 2\ / m 2 — u 2 

which solve ( p4"| ) and (f49|) as can be checked explicitly. It is not surprising at all that 
(|51~1) reproduces ( (29]) as cjf, — ► 0. 

Finally, in order to find gcr(x) we substitute ( |5*D| ) into (pE8"|). This is equivalent to 
solving 

[-d 2 x + m 2 - 2K 2 sech 2 [K(x - x )]} ^ (x) = (52) 

where ^oi x ) is defined in terms of gcr(x) in ([14]). For k 2 ^ m? the Schrodinger operator 
in ( |52"D has no normalisable zero-mode and supersymmetry is broken. Nevertheless, 
one can find non- normalisable solutions of the differential equation ([52]), which can 
be transformed into an hypergeometric equation [0], and extract gcr(x) in this way. 



The resulting ga(x) configurations are those found in |J. The specific gcr(x) given in 
corresponds to setting kx = f^Srzr i n (|50|) , namely, 



ga(x) = m + k I tanh 



KX 



1 111 + K 

-In 

4 771 — K 



tanh 



1 m + k 

kx + -In 

4 m — k 



(53) 



Note at this stage that we have yet to impose the saddle point condition fl22|). 
This will quantise uj^ and constrain it to the discrete set of values found by DHN Q 
as we now show. 

Substituting ( [46l) into the generic static saddle point condition ( ^2|) we obtain 

du 1 







(54) 



2tt (uj 2 — uj 2 ) \Jm 2 — u 2 

where we now integrate over uj along a contour C in the complex u plane, and the 
subindex n counts the number of fermions trapped in the single bound state of h_ 
produced by ga. The contour C is precisely the one used to define the Feynman 
propagator of a free Dirac particle of mass m, and in addition, it runs right below 
both simple poles of R- at uj = ±c<Jb||26||. 

The resulting static saddle point condition (|5"4"D is solved by requiring either that 



the term in the square brackets vanishes, or that the integral over frequencies vanishes. 
The differential equation resulting from the first possibility is 



ga" = 2ga{g 2 a 2 



m 



(55) 



which reproduces the kink configuration ([27]) discussed in the previous section with 
= 0. We thus focus on the other possibility, namely, that the integral in flo"4] ) 
vanishes. 

The integral in fl5"4]) is U.V. finite. We can therefore deform it by folding its uj > uj^ 
wing right on top of the remaining part of the real uj axis to the left of uj^. This can 
be further deformed into two circles wrapped around the simple poles at ±Ub and 
a "hair-pin" configuration wrapped around the left hand cut of picking up its 
discontinuity across the cut|26|. In other words, this integral picks up contributions 
from the completely filled Dirac sea (including the pole at uj = —uJb), where each 
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energy state (of the Dirac operator i$ — ga rather than hJ) is occupied by N fermion 
flavours, and the single bound state at u = +u b , which is occupied by n < N fermions. 

Because we have pulled out a factor of N in deriving @2]) we have to weigh the 
contribution of the simple pole of R_ at w = by n/N in ([54]) while the other 
contributions are weighed simply by 1 . Performing the integration, we find that 
the contribution from the cut is — m 2^ n2 g (l — ) while the poles at cu = ±u b yield, 
respectively, T m 2g in2e where, following @ we have defined 

u>b = m sinO and 



k = y m? — u% = m cos9. (56) 
Gathering all contributions to the integral (properly weighed), (151) yields 



duj 1 



c,n 2tt (ujI — lu 2 ) \Jm 2 — uj 2 m 2 sin26 



29 n 



namely, 



= % (57) 



which is precisely the result of |§. The simple poles in i?_ occur therefore at 

/ 7m 



u b = m sin i^—j . (58) 

One can now repeat the same analysis as in the previous subsection, of the par- 
tition function associated with fluctuations of fermions interacting with these static 
a(x) configurations. One finds again 0(2N) supermultiplets, where the n'th super- 
multiplet has mass M n = ^^sin (j^f) , (n < N) and contains all 0(2N) completely 
antisymmetric tensors of rank no < n, where Uq has the same parity of nH. 
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IV Conclusion 



In this paper we have solved the extremum condition on the effective action of the 
two dimensional Gross- Neveu model for static cr(x) configurations in the large N 
limit. Our method of calculation was direct, making use of elementary properties 
of one dimensional Schrodinger operators. The natural scale that appears in the 
Gross-Neveu model is that of its dynamically generated mass m. Because the latter 
is the natural scale of the CCGZ kink configurations as well, our derivation of these 
kink configurations was relatively simple and straight forward. It therefore may be 
considered as a clean and constructive proof that the CCGZ kinks are indeed static 
extrema of the effective action. We have also rederived the DHN extremal c{x) 
configurations in a very simple manner. To this end, however, we had to introduced 
their scale into the saddle point condition by hand, since it does not appear explicitly 
in the effective action. 

Our method may be applied also to a host of other two dimensional field theories, 
and in particular, to field theories that do not involve reflectionless static configura- 



tions, where inverse scattering methods are useless p7 
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been completed for their kind hospitality. I would also like to thank Joe Lykken for 
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Note Added in Proof: After submitting this paper for publication I have re- 



alised that Avan and de Vega |27| have already used the diagonal resolvent of a one 
dimensional Schrodinger operator to discuss solitons in "Large N" vector models. 
However, their discussion was limited to single particle quantum mechanics (0 + 1 
dimensional quantum field theory). Moreover, they did not make an explicit use of 
the Gel'fand-Dikii equation as was done here. I thank the referee for bringing this 
reference to my attention and for his useful remarks. 
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